Effective Viscosity of Confined Hydrocarbons 
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We present Molecular Dynamics (MD) friction calculations for confined hydrocarbon films with 
molecular lengths from 20 to f 400 carbon atoms. We find that the logarithm of the effective viscosity 
?7off for nanometer-thin films depends linearly on the logarithm of the shear rate: log7/ c ff = C— nlog-y, 
where n varies from f (solid-like friction) at very low temperatures to (Newtonian liquid) at very 
high temperatures, following an inverse sigmoidal curve. Only the shortest chain molecules melt, 
whereas the longer ones only show a softening in the studied temperature interval < T < 900 K. 
The results are important for the frictional properties of very thin (nanometer) films and to estimate 
their thermal durability. 
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The frictional and rheological properties of thin con- 
fined fluid films are of great importance in a number of 
engineering and scientific applications, e.g., in the con- 
text of lubrication. Thus, recent studies [1| have shown 
that for sliding contacts lubricated by organic or silicon 
oil bulk shear thinning can have a great influence on the 
friction at intermediate sliding velocities (mixed lubrica- 
tion). When wetting fluids are confined between solid 
walls at nanometer separation they often acquire solid- 
like properties and an increasing squeezing force is nec- 
essary in order to reduce the film thickness. If the solid 
walls are smooth the fluid molecules arrange in layers 
parallel to the solid walls and the squeeze-out occurs in a 
quantized way by removingone monolayer after another 
with increasing pressure [2tJ] ■ Sometimes the last one or 
two monolayers are so strongly bound that they cannot 
be removed by squeezing alone. Fluid films confined at 
the nanometer length scale exhibit viscosity enhancement 
and nonlinear flow properties characteristic of sheared su- 
percooled liquids approaching its glass transition [jj Q . 

Experiments for a large variety of fluids (including hy- 
drocarbon fluids and silicon oil) have shown that the 
logarithm of the effective viscosity r] c g for nanometer-thin 
films (typically 3 or 4 monolayers and contact pressures 
of the order of a few MPa) depends linearly on the log- 
arithm (here and elsewhere with 10 as the basis) of the 
shear rate: logoff = C — nlog7. If rj e g and 7 are mea- 
sured in SI units, for a large variety of fluids (at room 
temperature) C « 5 and n « 0.9. This linear relation 
has also been established in other experiments [H 0| and 
in computer simulations 1(J 11 1. 



We have performed a very extensive set of Molecular 
Dynamics (MD) simulations to probe the frictional prop- 
erties of thin layers of confined hydrocarbon molecules 
(with molecular lengths from 20 to 1400 carbon atoms). 
Some of these results may also be relevant for polymer- 
on-polymer systems as recent studies [13] have shown 
that in this case the shear deformations are localized to 
a band of material about 2.5 nm thick. Our results for n 
and C agree with the experimental observation at room 



temperature, but show that when the temperature in- 
creases n varies from 1 (solid-like friction) at very low 
temperatures to (Newtonian liquid) at very high tem- 
peratures. 

Our model is similar to those described in Refs. 
[lH, [13]. We consider a block and a substrate with atomi- 
cally flat surfaces separated by a polymer slab consisting 
of hydrocarbons with molecular lengths 20, 100 and 1400 
carbon atoms. The solid walls are treated as single layers 
of "atoms" bound to rigid flat surfaces by springs corre- 
sponding to the long-range elastic properties of 50 A thick 
solid slabs similar to our earlier papers (see, e.g. Refs. 
IB" The simulation box in x-y dimension is equal 
to 124.8 A x 124.8 A. In the following, periodic boundary 
conditions are assumed in the xy plane. 

Initially about half of the molecules are adsorbed on 
the block surface and half on the substrate surface. Two 
solids with adsorbed polymer slabs were put into contact 
and when the temperature was equal to the thermostat 
temperature everywhere we started to move the upper 
block surface. The temperature was varied from K to 
900 K to study the effect of temperature (and also melt- 
ing) on the shear stress. In our simulations, the polymer 
films are very thin (~ 3 nm), and the solid walls are 
connected to a thermostat at a short distance from the 
polymer slab. Under these circumstances we find that 
frictional heating effects are not important, and the ef- 
fective temperature in the polymer film is always close to 
the thermostat temperature. 

Linear alkanes C„H2,j+2 (with n = 20, 100 and 1400) 
were used as "lubricant" in the present calculations. The 
CH2/CH3 beads are treated in the united atom represen- 
tation [Ts|, [l9[ • The Lennard- Jones potential was used to 
model the interaction between beads of different chains 



U(r) = 46 



(i) 



and the same potential with modified parameters (ei, r\) 
was used for the interaction of each bead with the sub- 
strate and block atoms. 
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The parameters were eo = 5.12 meV for both the inte- 
rior and the end beads, and tq = 3.905 A. For the inter- 
actions within the C n H 2n +2 molecules we used the stan- 
dard optimized potentials for liquid simulations model 
HI 3' including flexible bonds, bond bending and tor- 
sion interaction, which results in bulk properties in good 
agreement with experimental data far below the boiling 
point 20]. Atomic mass 14 (for interior CH2 beads) 
and 15 (for the CH 3 end groups) were used. Within a 
C n H2n+2 chain we assume nearest neighbor C atoms are 
connected via springs with the spring constant k, which 
was chosen equal to 10 N/m. Note that this value is one 
order of magnitude smaller then the optimized 450 N/m 
[III], and was chosen such to facilitate a reasonable time 
step of 1 fs. We used an angle bending interaction of 
the form E(cos6)/kB — (l/2)/cbend(cos# — cos6*o) 2 with 
fc bcnd = 62543K and O = 2.0001 rad. For the dihe- 
dral interaction we used the functional form in term of a 
cosine Fourier series E(<j))/kB — Si=o c » cos l (</>) with pa- 
rameters c = 1009.99K, ci = 2018.95K, c 2 = 136.37K, 
C3 = — 3165.30K. Internal beads of separation greater 
than 3 units are treated similarly as beads from different 
chains. The number of molecules was equal to 1000, 200 
and 14 for the C20H42, C100H202 and C1400H2802 systems 
respectively. The hydrocarbon films at room tempera- 
ture consisted of 6 to 8 monolayers of molecules between 
the solid surfaces. The (nominal) squeezing pressure po 
was usually 10 MPa. 

We have chosen the polymer-wall atom bond to be so 
strong that no slip occurs at these interfaces. This is 
the case with n = 2.92 A, t\ = 160 meV. The lattice 
spacings of the block and of the substrate are a = b = 
2.6 A. 

If v is the sliding velocity and d the film thickness, we 
define the shear rate 7 = v/d and the effective viscosity 
r? e ff = crf/7, where at is the frictional shear stress. Fig. Q] 
shows the logarithm of the effective viscosity as a function 
of the temperature [21| for the (a) C20H42, (b) C100H202 
and (c) C1400H2802 system, at four sliding velocities (from 
top to bottom) 0.3, 3, 10 and 100 m/s. Increasing the 
velocity results in a reduced effective viscosity, i.e., the 
thin films exhibit shear thinning. 

The velocity dependence of the MD data in Fig. Q]can 
be very well fitted by the formula rj c tt = Bj^ 11 or 



logoff = C - nlogj, 



(2) 



where C = logB. This is illustrated in Fig. [2] for two 
cases. The rest of the cases also show that a linear fit 
is appropriate with the exceptions of the results at very 
high temperatures. The C100H202 system shows the high- 
est effective viscosity due to entanglement. The C20H42 
is nearly liquid above 300 K and shows less entanglement 
whereas the sliding in C1400H2802 takes place at nearly 
one interface (see figure [6]). We find that the parameters 
C and n depend on the temperature. The data points in 
Fig. |3] show the temperature dependence of the index n, 
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FIG. 1: (Color online) The logarithm of the effective viscosity 
as a function of the temperature for (a) C20H42, (b) C100H202 
and (c) C1400H2802 system at four sliding velocities: 0.3 m/s 
(■), 3 m/s (A), 10 m/s (♦) and 100 m/s (o). 
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FIG. 2: (Color online) The data points for the logarithm 
of the effective viscosity as a function of the logarithm of the 
shear rate at four sliding velocities: 0.3, 3, 10 and 100 m/s. 
The lines represent the linear fits to these data points. The 
systems are C100H202 at 600 K and C1400H2802 at 300 K. 
The values of the slope and the intercept of the C100H202 
are -0.520 ± 0.038 and 2.260 ± 0.351 and for the C1400H2802 
system -0.958 ± 0.024 and 6.753 ± 0.228. 



while the solid lines are fits to the data points using the 
inverse sigmoidal curve: 



1 



1 + (T/T c ) a ' 



(3) 



The parameters a and T c are given in table HI At low 
temperatures n — 1 as expected for dry friction. That 
is, the frictional shear stress at = ?7cff7 = 10 7 1_ ™ is 
independent of the shear rate when n = 1 as expected for 
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FIG. 3: (Color online) Values of the parameter n in equation 
(2) as a function of the temperature for the systems C20H42, 
C100H202 and Ci4ooH2802- The solid lines are fits to the nu- 
merical data using equation (3) with the parameters a and 
T c given in table [I]. 





C20H42 


C100H202 


C1400H2802 


T c (K) 


353 ± 28 


642 ± 22 


840 ± 30 


a (-) 


4.09 ± 1.24 


3.68 ± 0.46 


2.79 ± 0.30 



TABLE I: Table showing the T c and a in equation (3) for 
each system. The standard deviation of the fitting parameters 
is also indicated. 
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FIG. 4: (Color online) The C parameter in equation (2) 
also follows a sigmoidal curve when the temperature is varied. 
This figure shows the dependence of C on n. All three systems 
C20H42, C100H202 and C1400H2802 are represented. 




dry friction at low temperatures (no thermally activated 
creep). At high temperature n approaches as expected 
for a Newtonian fluid where the frictional shear stress is 
proportional to the shear rate. 

The parameter C depends on the units used for r] e s 
and 7 and here we assume SI units. Remarkably we find 
a linear relation between C and n for all the systems and 
temperatures we have studied, see Fig. H) The extrapo- 
lation of C to n = gives C « —3.8 ± 0.2. This gives the 
effective viscosity r] e e = _L0 C = (1.6 — 2.5) x 10~ 4 Pas. Ex- 



periments have shown [22j that for a wide range of fluids 
the fluid viscosity at the boiling point is « 2.2 x 10~ 4 Pas, 
i.e., nearly the same as we deduce for confined fluids when 
extrapolating our data to n = (Newtonian fluid). This 
result is extremely interesting but perhaps not entirely 
unexpected, as n — corresponds to high temperatures 
where the separation between the solid walls (see Fig. [5] 
in which the sliding velocity is 10 m/s) is much larger 
than for lower temperatures and where the mobility of 
the fluid molecules may be similar to that in the bulk 
fluid close to the boiling point. The C100H202 film has 
a smaller film thickness than the other systems at low 
temperatures due to the fact that it only has 6 monolay- 
ers at these conditions. The other two systems have 7 
monolayers at the same conditions and have thus thicker 
films. This matter has been discussed in reference 14 j. 

Fig. |3] shows that the longer the hydrocarbon chain 
length the higher the temperature necessary for the film 




300 600 
temperature (K) 



900 



FIG. 5: (Color online) The polymer film thickness as a func- 
tion of the temperature. The sliding speed is 10 m/s. All three 
systems C20H42, C100H202 and C1400H2802 are represented. 



to behave as a Newtonian fluid. This is also illustrated 
in Fig. |6] which shows the relative number of atoms trav- 
eling at velocities below v, as a function of v, for all three 
polymer films at 600 K. It can be seen that the velocity 
gradient is going from nearly Newtonian in the case of 
C20H42 to solid-like in the C1400H2802 case. 

Note that for the hydrocarbons C100H202 and 
C1400H2802 at room temperature n ss 0.9, as also found 
experimentally for many confined fluids. When n « 0.9 
from Fig. @] we get C ~ 6 which is a little larger than 
found experimentally (at room temperature) C exp ~ 5. 
This may reflect somewhat different confinement condi- 
tion, e.g., differences in contact pressures (10 MPa in our 
study as compared to a few MPa in most of the experi- 
mental studies). However, in both cases the contact pres- 
sure is so small that one does not expect any significant 
dependency of the shear stress on the contact pressure 
[l3j, except if there is a pressure induced change in the 
number of confined layers. 

It is well known that the viscosity of fluids at high pres- 
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FIG. 6: (Color online) The relative number of atoms moving 
at the velocity below v, as a function of v, at 600 K. It can be 
seen that the velocity gradient is going from nearly Newtonian 
in the case of C20H42 to solid-like in the C1400H2802 case. 



sures may be many orders of magnitude larger than at 
low pressures. Using the theory of activated processes, 
and assuming that a local molecular rearrangement in a 
fluid results in a local volume expansion, one expects an 
exponential dependence of viscosity -q on the hydrostatic 
pressure p, r\ = r]oenp(p / po) , where typically (for hydro- 
carbons or polymer fluids) po f» 10 8 Pa (see, e.g., Refs. 
Hill). Here we are interested in (wetting) fluids con- 
fined between the surfaces of elastically soft solids, e.g., 
rubber. In this case the pressure at the interface is usu- 
ally at most of the order of the Young's modulus, which 
(for rubber) is less than 10 7 Pa. Thus, in most cases 
involving elastically soft materials, the viscosity can be 
considered as independent of the local pressure. 

One of us have recently studied rubber friction on 
rough surfaces |25(. For unfilled styrene butadiene rub- 
ber we found the transfer of a thin smear film to the 
substrate. In this case the shear during sliding may be 
localized to a thin (a few nanometer) interfacial layer 
which may exhibit frictional properties very similar to 
what we have observed in our simulations. Indeed the ex- 
perimental data indicated a frictional shear stress of the 
form predicted above with n « 0.91 and 10 c w 1.3 x 10 5 
(in SI units), in close agreement with the result of our 
simulations. 

To summarize, we have presented results of molecular 
dynamics calculations of friction performed for a block 
sliding on a substrate separated by k 3 nm thick poly- 
mer films where the alkanes had 20, 100 and 1400 carbon 
atoms. In all cases we found that the logarithm of the 
effective viscosity is proportional to the logarithm of the 
shear rate, logoff ~ C — nlog7. The index n varies from 
1 (solid-like friction) at very low temperatures to (New- 
tonian liquid) at very high temperatures, following an in- 
verse sigmoidal curve. The C parameter is proportional 
to n and as n — > 0, 10 c extrapolates to the viscosity of 
the bulk fluid at the boiling point. At room temperature 



the parameters n and C have been found to be close to 
what has been observed experimentally for a large num- 
ber of fluids. 
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